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Abstract 

We study asymptotic expansions of spin-spin correlation functions for 
the XXZ Heisenberg chain in the critical regime. We use the fact that 
the long-distance effects can be described by the Gaussian conformal 
field theory. Comparing exact results for form factors in the XYZ and 
sine-Gordon models, we determine correlation amplitudes for the leading 
and main sub-leading terms in the asymptotic expansions of spin-spin 
correlation functions. We also study the isotropic (XXX) limit of these 
expansions. 



1 Introduction 



In the domain of two-dimensional integrable models, it is, in general, still a challenging prob- 
lem to compute correlation functions in the form of compact and manageable expressions. 
For lattice systems, a few methods of computation have been developed f^-Q]: in particular, 
it is possible in some cases to obtain exact integral representations of correlation functions. 
However it is still quite difficult to analyse those expressions, and especially to extract their 
long-distance asymptotic behavior. 

On the other hand, at the critical point, when the gap in the spectrum of the lattice 
Hamiltonian Hiatt vanishes, the correlation length becomes infinite in units of the lattice 
spacing. As a result, the leading long-distance effects can be described by Conformal Field 
Theory (CFT) with the Hamiltonian Hcft fH- In general, it is expected that the critical 
lattice Hamiltonian admits an asymptotic power series expansion involving an infinite set of 
local scaling fields (see e.g. P, 

Hiatt ~ Hcft 



+ Ai j dxOi + \2 e'^^'^ j dx02 + ... . (1.1) 



Here dk denotes the scaling dimension of the field Ok, and the explicit dependence on the 
lattice spacing e is used to show the relative smallness of various terms. All fields occurring 
in (^]^) are irrelevant and the corresponding exponents of e are positive. Such an expansion 
also involves a set of non-universal coupling constants depending on the microscopic 
properties of the model, and whose values rely on the chosen normalization for Ok- 

The asymptotic expansion (|1.1| ) is a powerful tool to study the low energy spectrum 
of lattice theories. At the same time, in order to analyse lattice correlation functions, one 
should also know how local lattice operators Oiatt are represented in terms of scaling fields. 
Just as the Hamiltonian density, they can be expressed as formal series in powers of the 
lattice spacing e, 

Olatt ~ Cm, e'^-^ Om, + Cm, s'^"'^ Om, + ■ ■ ■ , (1-2) 

where again the non-universal constants Cm depend on the normalization of the scaling 
fields. 

The knowledge of the explicit form of ( |1.1D and ( |1.2D enables one to obtain the long- 
distance asymptotic expansion of lattice correlation functions in powers of lattice distances. 
For example, in the case of a vacuum correlator, the leading term is given by the CFT 
correlation function of the operator occurring in (|1.2| ) at the lowest order in e, whereas 
subleading asymptotics come from both higher order terms in ( |1.2D and from the perturbative 
corrections to the CFT vacuum toward the lattice ground state. In this expansion, the exact 
values of the exponents follow from the scaling dimensions of the CFT fields, while the 



corresponding amplitudes come from the values of the constants appearing in (IT) and 
([1.21) . Thus, to investigate quantitatively the long-distance behavior of lattice correlation 
functions, one needs to determine (for a fixed suitable normalization of the scaling fields) 
exact analytical expressions for the non-universal constants A^ and Cm- 

The aim of this article is to study expansions of the type ( |1.2p in the case of the XXZ 
spin 1/2 Heisenberg chain. In particular, we calculate the first constants Cm occurring in 
the expansion (|0) of local spin operators. 
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The article is organized as follows. 

In Section ^, we recall the definition of the XXZ spin chain (see e.g. |P for details), and 
discuss its continuous limit, the Gaussian CFT P-p!l|. We review in particular how, from 
the analysis of the global symmetries of the model, one can derive selection rules for the 
expansions (1.2) of the lattice operators in terms of the scaling fields of the Gaussian model. 
This enables one to predict the structure of the long-distance asymptotic expansions of the 
spin-spin correlation functions. 

The problem of computing the constants occurring in the expansions ( |1.2D of spin opera- 
tors is the subject of Section ^. There we explain how, moving slightly away from criticality, 
and comparing exact results obtained for the XYZ and sine-Gordon models, one can quanti- 
tatively connect lattice spin operators to scaling fields. This gives us access to the correlation 
amplitudes of the spin-spin correlation functions. Our predictions are gathered in Section ^, 
where they are compared to existing numerical data. 

In Section ^, we study these spin-spin correlation functions in the isotropic (XXX) limit 
by means of the exact Renormalization Group (RG) approach. This section contains an 
erratum of Section 5 of |12|. 

Finally, we conclude this article with several remarks. 



2 XXZ spin chain in the continuous Hmit 



2.1 Preliminaries [§,|T^ 



To illustrate the problem of the determination of correlation amplitudes, we consider an 
example, the XXZ Heisenberg chain of spins 1/2. It is defined by the Hamiltonian 

J oo 

Hxxz = - 2 E {'^r^r+i + + ^^f^f+i} ' (2-1) 

l=—oo 

where the spin operators (a = x, y, z) denote the conventional Pauli matrices associated 
with the l-ih site of the infinite lattice, and A is an anisotropy parameter. As Hxxz(— 
can be obtained from Hxxz (<^) A) by a unitary transformation, we choose in the following 
the coupling constant J to be positive. 

The nature of the spectrum of the infinite chain depends on the value of the anisotropy 
parameter A. The study we present here concerns the critical regime of the chain, which 
corresponds to the domain — 1 ^ A < 1. We shall use the parameterizations 

A = cos(7rr?) (0 < < 1) , (2.2) 

and 

J = ^.""^ . (e > 0) . (2.3) 

esm(7rr?) 

To describe leading long-distance (low-energy) effects, it is useful to consider a continuous 
limit of the lattice model. This limit can be obtained from the representation of ( p.l| ) in 
terms of lattice fermions, 



Hxxz = -J E{ WV'i+i + 4+1^1 + y (1 - 2^^/) (1 - 2W+iV'^+i) } , (2.4) 
I 
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where the fermionic operators are related to the spin operators through the Jordan- Wigner 
transformation, 



af = 1-2^^,^/, W = n^r^r' ^i = ll''r''t^ (2-5) 

and = {af ziziaf)/2 . With the parameterization ( |2.3D of the constant J, the spin-1/2 spin 
wave dispersion relation has the form [13|, 14] 



£{k) = -'^. (2.6) 

The ground state of the chain has all levels filled with \k\ < 7r/2. Linearizing the dispersion 
relation in the vicinity of the Fermi points kp = ±7r/2, one can take the continuous limit 
of the model in which the lattice operators ipi are replaced by two fields, iPr{x) and tpiix), 
varying slowly on the lattice scale: 

ipi (X i)R{x) + e~''^ ipLix) . (2.7) 

Here x = le, and thus the parameter e can be interpreted as a lattice spacing. The continuous 
Fermi fields are governed by the Hamiltonian 



H 



Thirring 



2tt 



A precise relation between the anisotropy parameter A and the four-fermion coupling con- 
stant g depends on the choice of the regularization procedure for the continuous Hamiltonian 



(2.8) and is not essential for our purposes. What is important is that the quantum field the- 



ory model ( |2.8| ) (which is known as the Thirring model |15, 161) is conformally invariant and 



equivalent to the Gaussian CFT (see the preprint collection for a historical review of 
bosonization). 

2.2 Exponential fields in the Gaussian model 

The Gaussian model is defined in terms of one scalar field if, which satisfies the D'Alembert 
equation^] 

{df -v'^d'^)ip = 0, with?; = l, (2.9) 
and the boundary condition 

dt^pix,t)\x^±oc = dx^pix,t)\x^±oo = . (2.10) 
Assuming that the equal-time canonical commutation relations are imposed on the field, 

[(/?(x), dt(p{x') ] = Sni d{x — x') , 



^ The value of the spin-wave velocity v is determined from the slope of the dispersion relation p.q ) at the 
Fermi points and from the identification of x with I e. 
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one can write the Hamiltonian of the model in the form 

Hcau! 

where 



y ^{T/j + Tl} + const, (2.11) 



Tr(x) = ^{d,v - dtvf , Tl{x) = ^{d,^ + dtvf ■ (2.12) 

lb lb 

As is usual in CFT , it is convenient to set a class of conformal primary fields among 
all scaling fields. In the case of the Gaussian model, these conformal primaries include right 
and left currents, 

id,-dt)^, {d,+dt)ip, (2.13) 

along with exponential field^ |l^. With the boundary condition ( 2.1(]| ), the latter can be 
defined as 

Os,nix, t) = A'^'-" Ihn^ exp | ^ j dx d^x' , t) | exp | ip{x + e, t) | , 

(2.14) 

where s, n are integers and 

2 2 

s r] n 

dsn = — - + —■ 2.15 

The regularization parameter A, which has the dimension [length]"^, is introduced in the 
definition ( ^.14 ) in order to provide a multiplicative renormalization of the fields. Notice 
that Os^n obey the simple Hermiticity relation 

Oj „ = . (2.16) 

To completely define these exponential fields, we should also specify some condition which 
fixes their multiplicative normalization. By a proper choice of A in ( 2.14| ), one can impose 
the following form of the causal Green's functions in the Euclidean domain x"^ — > 0: 

, sn 

{ T a,n(x, t) Oj „(0, 0) ) = (^) ' ix' - t^)-<- . (2.17) 



We will later refer to Eq. (2.17) as the "CFT normalization condition". 
2.3 Global symmetries 

To draw a precise link between the Gaussian CFT and the XXZ spin chain, it is important 
to examine and identify their global symmetries. 

Let us first consider the Gaussian model. The Hamiltonian ( 2.11 ) is manifestly invariant 
under the U{1) rotations, 

VaifV-^ = ip + ^ , U„a,nU-i = e^^°a,n , (2.18) 



^Here we are not considering the orbifold (Ashkin- Teller) sector of the Gaussian model. 
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where the operator Uq, can be written in the form 



oo 



The CFT model is also invariant under the parity transformation P : ip{x, t) — > (p{—x, t), the 
time reversal T : Lp{x,t) (p{x,—t), and the reflection C : (p ^ —ip. Using the definition 
( 2.14| ), one can show that these transformations act on the exponential fields as 



P Os,nix, t)F = e-— U^n a,-n(-x, t) , (2.20) 

TOs,n{x,t)T = 0-s,n{x,-t), (2.21) 

Ca,nC = 0_,,„„(x,t). (2.22) 

Note that the reflection C can be naturally considered as a charge conjugation in the theory. 
As usual, C and P are intrinsic automorphisms of the operator algebra, contrary to the 
anti-unitary transformation T which acts on c- numbers as follows, 

T (c-number) T = (c-number)* . 



Let us now identify the above transformations with the global symmetries of the spin 
chain. First of all, it is natural to define the action of the U{1) rotation on the lattice as 
follows, 



U„ af = e±*" af , U« af U^' = af , (2.23) 

where 

with S, = IY.''j- (2.24) 



U„ = e 



ictSz 



1 



Such an identification has an important consequence. Indeed, as the two rotations of angle 
27r and —2ti are indistinguishable trivial transformations in the lattice theory, we should set 

U2,r = U_27r = i^I with u = ±1 (2.25) 



in the corresponding Gaussian model. The condition ( 2.25| ), along with ( |2.18| ), implies that 
ip has to be treated as a compactified field: 

47r 

(p = LP ^ — Z . 



According to ( |2.24| ), the sign factor u = +1 in Eq. (2.25) occurs for states with an integer 



eigenvalue of the operator Sz- The corresponding linear subspace of the whole Hilbert space 
can be constructed from the thermodynamic limit of finite chains with an even number of 
sites, and will therefore be referred to as the "even sector". Similarly, the condition u = —1 
defines another linear subspace which will be called the "odd sector". 

The actions of charge conjugation and time reversal can be naturally identified, in the 
lattice theory, with the following transformations, 

CafC = af , C (jf C = -af , (2.26) 

TafT = af , T T = -af . (2.27) 
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Let us recall here that the time reversal is an anti-unitary transformation, so that C and T 
correspond to different symmetries of the XXZ chain, even thought they act identically on 
the spin operators. 

As for the parity transformation, its action in the lattice model depends on the choice of 
the sector specified by the sign factor v in ( |2.25 ). Indeed, v = +1 implies that the considered 



infinite chain is defined as the thermodynamic limit of finite lattices with an even number of 
sites. Such finite lattices clearly do not possess any invariant site with respect to the parity 
transformation. Therefore, 

Yo1'^ = ol_i for zv = +l, (2.28) 

whereas the "naive" action of P is valid only in the odd sector: 

PafP = cj^i for v = -\. (2.29) 

2.4 Selection rules 

The global symmetries and the knowledge of their action on lattice and continuous operators 
provide selection rules for the set of scaling fields which can occur in expansions (0) and 
(|l.2D . Here we examine what are these selection rules for the conformal primary fields in the 
expansions of the spin operators. 

First, since the whole operator content of the Gaussian model is given by the primary 
fields described above and by their conformal descendants, it is easy to see that O±i,o are 
the local fields with the lowest scaling dimension which may occur in the expansions of cr^ : 

4~C7oe''i'0 O±i,o(0) + ... . (2.30) 

Furthermore, assuming that the conjugation in the lattice theory is defined in such a way 
that o"" (a = x,y,z) are Hermitian operators, we conclude from ( 2.16 ) that the constant Cq 
in ( p.30| ) is real. 

Let us now consider the expansion of the lattice operator ctq. Due to the ^7(1) invariance, 
it can contain only the primary fields {dx ±dt)(p and Oo^n, n £ Z, along with their conformal 
descendants. Using definition ( |2.14| ), it is easy to show that the fields Oo,n with odd n are not 
mutually local with respect to O±i,o- As the latter are the leading terms in the expansions 
(|2.30| ) of a^, it would contradict the mutual locality of the lattice operators and af if 
there were any 00,2^+1; m £ Z, in the series for Uq. From the C, T invariances and the 
Hermiticity of ctq, one can moreover predict that the primary fields are allowed to appear 
only as linear combinations of dtip and i (Oo,2m — Oo,-2m) with real coefficients. Finally, let 
us consider the parity transformation. According to equations ( 2.20| ) and ( |2.25| ), 



P a,2m(x, t)F = l^"' Os,-2m{-X, t) , (2.31) 

and thus one has to examine even and odd sectors separately. In the odd sector = — 1), 
the parity (|2.29| ) prohibits the presence of the fields Oo,4m — Oo,-4m with m G Z. Therefore 
the expansion has to be of the form. 



ag-eCo^a^O) 
^ +00 

+ ^ (^o,4m-2 - Oo,-4m+2)(0) + descendants . (2.32) 
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But, since ( p.32| ) is a local operator expansion, it cannot depend on the choice of the sector 
in the Hilbert space. Thus, ( 2.32| ) should be valid in the even sector as well. Owing to this 
and to Eqs. ( 2.28| ), ( p. 31 ), we can determine the action of the lattice translation, 

KafK-^ = af^,, (2.33) 

on the primary fields 00,401-2: 

K Oo,4m-2ix) = -Oo,4m-2(X + e) . (2.34) 

With this relation, the expansion ( |2.32| ) can be written in the more general form, 
2tt^ 



+ Yl (^o,4m-2 - Oo,-4m+2) (x) + descendants , (2.35) 



m=l 



where x = le. Notice that the exact value of the coefficient Cg in ( 2.35| ) is determined from 
the comparison of Eqs. ( 2.ig| ) and ( |2.24| ). 

Using the same line of arguments, one can extend the expansion ( |2.30| ), 

-. +00 

^f-l^Y. (-1)'" (^±1,2- + 0±i,-2m) {x) + descendants , (2.36) 

m=0 

and determine the action of the lattice translation K on 0±i^2m'- 

K 0±i,2m(x) = (-1)™ 0±i,2™(x + e) . (2.37) 

Since all exponential fields 0±s^2m{x) with integers s and m can be obtained by means of 
operator product expansions of the fields 0±i^2m{x) and 00,401-2(2^)5 we deduce from ( |2.34|) 
and ( |2^ that 

K 0,,2™(x) = (-1)™ 0,,2™(x + e) . (2.38) 

Notice that, in the process of bosonization, the fermionic fields tpji and ipL ( |2.7|) ar e identified 
with the exponential fields Oi,_i and Oi,i respectively. Hence equation ( p.38|) is indeed 
consistent with 



It is not difficult to extend our symmetry analysis to derive the selection rules for expan- 
sion ( |L^ ) of any local lattice operator. This procedure can, in particular, be applied to the 
case of the Hamiltonian density. As a result, the following form of the low energy effective 
Hamiltonian for the XXZ spin chain is suggested: 



Hxxz ~ Hcauss + ^\ Yl ^me'^"'"'""^ (Oo,4m + Oo,_4m) (a:^) + descendants \ . 

^ m=l 

(2.39) 

As was mentioned in Introduction, it is necessary to choose the normalization of the 
scaling fields to give a precise meaning to the couplings in ( 2.39 ) as well as to the real 
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constants C^, in expansions ( p. 35 ) and (2.36). In this paper we adapt the CFT normal- 
ization ( 2.17 ). With this normahzation condition, the first couphng constant Ai in ( |2.39| ) 
has been obtained in |jl^, together with the leading contributions of descendant fields (see 
also Ref. for a qualitative analysis of the descendent field contributions): 



H 



xxz 



H 



Gauss 



(T| + T 



+ 



where 



Ai 

A+ 
A_ 



4r(i/7?) 

r(i - l/r?) 



r 1 + 



2~2riJ 



20Fr(i + 



2/r,~2 



— tan , 
27r V2 



2r] 



1] ^{2^)^^i2^) 

i2vrr(^)r3( 1 ) 



(2.40) 

(2.41) 
(2.42) 
(2.43) 



2.5 Vacuum spin-spin correlation functions 



At this stage, without the knowledge of the constants occurring in expansions ( p.35 ) and 
( 2.36| ), it is already possible to predict^ the exact values of the exponents for the vacuum 
spin-spin correlation functions. They follow immediately from the scaling dimensions of the 
fields occurring in (|2|3|), ( p6| ) and 1^^). 

It may be worth recalling at this point, that the vacuum sector of the infinite XXZ chain is 
infinitely degenerate. In general, the boundary conditions imposed on the finite critical chain 
do not preserve all the global symmetries discussed above, and they may be spontaneously 
broken at the thermodynamic limit. Here we consider only translational invariant vacuums 
with unbroken parity: 



vac , 



vac , 



vac , 



(2.44) 



To fulfill these requirements, we shall treat the infinite XXZ chain as the thermodynamic 
limit of finite chains subject to periodic boundary conditions. Then, different vacuum states 
can be distinguished by means of the operator Sz (|2.24|) : 

S'^|s) = s|s), 2s £Z. (2.45) 

Also, by a proper choice of phase factors, one can always set up the conditions 

C|s) = |-s) and T|s) = |-s). (2.46) 

In the continuous limit, the vacuum | s) flows toward the conformal primary state with right 
and left conformal dimensions equal to dsfi/2, where ds^n is given by (|]l|) loHll. This 



^It has been done for the first time in Ref. j^. 
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implies in particular that, for the XXZ chain with a finite number of sites ^ 1, the 
difference of vacuum energies corresponding to the states | s ) and | ) is 2-Kdsfi/N + o{N~^). 

From our previous analysis, one can now predict the following asymptotic expansions for 
the time-ordered correlation functions: 



(raf+/t)cTj(0)) 



A 



(-i)'i ri//+ 



'7 I 1 



+ 



1 



1 



B 



+ 



B 



(l+l-) 



— + o(r2 log M'-'/^)} + ... 



+ 



(-1)' ^ 



TT^r] l+l^ t 2 V /_ / 
2-77 / /+ Z_ 

1 + ^ ( 7^ + 



4(l-r?)VL Z+ 
5 



+ 0{l-^logl,l''-'/^) 
1? J 



where 



(2.47) 



(2.48) 



Z± = / ± - > 1 



(2.49) 



The coefficients in the asymptotic expansions ( 2.47 ) and ( ^.48 ) do not really depend on the 
choice of the vacuum state \s), thus any finite s G Z/2 can be chosen for the averaging 

/ V ^ {s\...\s) 
\---/ - {s\s) ■ 

In Eqs. ( p.47| ), ( p. 48 ), the correlation amplitudes A, A and A^ are simply related to the 
first constants occurring in the expansions ( p.35 ) and ( |2.36 ): 



^ = 2(Co)^ i = (Ci)^ 



A 



I? 



(2.50) 



They will be computed in the next section. At the same time, the constants B, B, Bz and 
Bz appearing in ( pTl) and (pip can be determined by methods of conformal perturbation 
theory |23] based on the effective Hamiltonian ( 2.40| ). In particular, the constant B was 
obtained in Ref. ll^ from second order perturbative calculations. 



B = ^ 



27r^ 



16 i sin2(27r/r/) (1 -r/)(2 -r/) 



V^'(l/r/)-V'(3/2-l/r?) , 



(2.51) 



where ip'{^) = d'^logT(z) and Ai is given by ( 2.41 ). The constant B^ can be computed 
similarly: 



Bz = ^ 



V 



4 (r/-2)2 ■ 



(2.52) 
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On the contrary, the determination of B and i?^ do not require calculations beyond the first 
order. They read explicitly, 

B = {1- rif . (2.54) 

3 Calculation of correlation amplitudes 

The purpose of this section is to explain how the local operators af and erf can be quanti- 
tatively related to the scaling fields ( 2.13| ) and ( p. 14 ), that is, how one can compute explicit 



analytic expressions for the constants Cm and occurring in ( 2.3E| ) and ( |2.36| ) for the fixed 



CFT normalization ( p.l7| ). We concentrate here on the first terms of the expansions ( ^.35| ), 



(2.36), which provide the main asymptotic behavior of the spin-spin correlation functions 



(P^ ) and (|1|): 



af ~ Co £2 0±ifl{x) + (-1)' ^ e^+^ (0±i,2 + ©±1,-2) {x) + . . . , (3.1) 
dMx) + (-1)' ^ (Oo,2 - 00,-2) (x) + . . . . (3.2) 



27r^ ' ' ' 2i 



Note that in ( |3.2| ) each of the two terms is either leading or sub-leading according to the 
value of r/, i.e. of the anisotropy parameter A, whereas in ( p.lj ) the second term is always 
subleading. Nevertheless, this term gives rise to noticeable numerical corrections to the lead- 
ing asymptotic behavior of correlation functions and has already been studied numerically 
(see Section Q) . 

3.1 Even sector of the infinite XYZ chain 

To find the quantitative relation between the spin operators and the local scaling fields, it 
is useful to move slightly away from criticality and, instead of ( |2.1| ), to consider the XYZ 
chain, 

^ 00 

HxYZ = "2 E ^r+i + Jy ^f<i + Jz ^f+i} • (3.3) 

Z=— 00 

Without loss of generality we assume here that J^. > Jy > \Jz\- The XYZ deformation has 



a remarkable feature: it preserves the integrability of the original theory [|,24]. Nowadays, 
the structure of the Hilbert space of the infinite XYZ chain is well understood. We recall 
here some basic facts that will be useful for our analysis. 

In the case of the XYZ spin chain, the global symmetry group discussed in the previous 
section is explicitly broken to the subgroup generated by the lattice translation K, the C, P, 



T transformations and the rotation ( 2.23| ) with q = vr: U^. In this section, we concentrate on 



the even sector of the XYZ spin chain (defined as the thermodynamic limit of finite chains 
with an even number of sites), which implies the condition 

Vl = l. (3.4) 
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Let us denote by Tig the eigenspace of the operator Sz corresponding to a given eigenvalue 
s G Z. Then, the Hamiltonian ( |3.3D acts as, 

HxYZ : T~(-s-2 © © , 

and therefore the infinite degeneracy in the vacuum sector of the XXZ chain is reduced to 
the two states 

I S )XYZ G ©r=-oo W^+2fc ( S = 0, 1 ) 

satisfying the condition^ 

Un\s)xYZ = e'^' \S)XYZ ■ (3.5) 

This requirement, along with the conventional normalization of vacuum states, (vac | vac ) = 
1, defines | s )xYZ up to an overall complex phase factor. Since the time reversal transforma- 
tion acts on states as the complex conjugation, one can eliminate (up to sign) the ambiguity 
of such a phase by imposing the condition 

T I S )xYZ = I S )xYZ . 

To gain physical intuition about the vacuum sector, it is useful to consider a limiting case 
where the Hamiltonian ( |3.3| ) simplifies drastically. For Jy = = 0, the vacuum sector of the 
finite periodic XYZ chain contains two pure ferromagnetic states | vac (j = 0, 1). With a 
proper choice of the overall phases of these states, one may always set up the conditions 

I vac ) = I vac ) ^^^^'^ , T | vac ) (-^'^ = | vac ) ^^'^ . (3.6) 



Moreover, since the charge conjugation matrix ( 2.26 ) can be identify with H^crf, one has 

C I vac = I vac) (^'^ . (3.7) 

When the couplings Jy and Jz are non-vanishing, the pure ferromagnetic states are no longer 
stationary states, but it is still possible to introduce two vacuums |vac)^-'^ satisfying 



and (3/7) through the relation 

|s)xYZ = ^ {|vac)W + (-l)^|vac)W} (^ = 0,1) . (3.8) 

The Hilbert space of the XYZ chain contains two linear subspaces V^-'^ (j = 0, 1) as- 
sociated with the vacuums \vac)^^h In the spectrum of the model, there exist kink-like 
"massive" excitations, B+ and B_, such that the corresponding Zamolodchikov-Faddeev 
operators intertwine these subspaces V^-'^: 

One can therefore generate two sets of asymptotic states in the form 

B,,„(A;2^)...B,,(fci)|vac)(^') G V(^') (m = 0, 1, 2 . . . ) , (3.9) 



* The corresponding energies for a chain with a finite number of sites N ^ 1 are asymptotically 

.(1) _ p(0) 



degenerate in the sense that E)^' - E)^' = 0(e-'=°"=*™) (see e.g. Ref. H). 
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where the Zamolodchikov-Faddeev operators Bo-(A;) (a = it) depend on a quasi-momentum 
k: 

KB±{k)K-^ = e''' B±{k) , [UxYz,B±{k)] =£{k)B±{k) . (3.10) 

The dispersion relation £ = £{k) of the fundamental excitations was calculated in work p^ . 



It is recalled in Appendix A (see Eqs. ( A.8 ), ( |A.9| )). The operators B-t satisfy also the 
conditions: 

c Bi^-^'^') c = T {-ly b5!-^'^') , B^!"^'^') = ± (-i^ b^'^-^'^ . (3.11) 



In Eqs. ( 3.11 ), B^ ' denotes the restriction of the operator B± when it acts on the 
subspace V^-'^. 

Any local lattice operator Ciatt (i-e- any operator which can be written as a local com- 
bination of af) leaves the subspaces V^-'^ invariant: 

Olatt : V(^) ^ V(^) . 

Furthermore, the algebra of local operators ^loc acts invariantly on the component of the 
subspace V^-^^ generated by the states (|3.9D . For < 0, the kinks B± do not produce bound 
states, and the two sets of asymptotic states (p.9[) obeying the condition 



-7r/2 < A;i < ^2 • • • < ^2m < 7r/2 

form complete in-bases in these unitary equivalent spaces of representation of A\oc- At first 
glance, to construct explicitly the representations of .4ioc, one should put at one's disposal 
the whole collection of in-basis matrix elements for an arbitrary local operator Ciatt S ■A.\oc- 
As a matter of fact, using the so-called crossing symmetry (see e.g. [|2|,^), one can express 
all possible matrix elements in terms of those of the form 

(^■)( vac I Oiatt B,,„ {k2^) . . . B,, (fci) I vac )(^') = Outt I B,, {k^) . . . B,,^ (^2™) . 

(3.12) 

Such matrix elements are known as form factors. 

Currently, there exists a formal procedure which allows one to express form factors of 
local operators in terms of multiple integrals [26,27]. Unfortunately, it is difficult to apply 
in practice, even in the case of the local spin operators themselves and for form factors 
involving only a small number of excitations. For our purposes, we merely need the explicit 
form of Vacuum Expectation Values (VEV) and of two-particle form factors. From the 
relations ( ^.44 ), (|3.6|) and (^), it follows immediately that the VEVs of and vanish, 
whereas 

i^){at) = {-iy F , (3.13) 

where F depends on the two ratios Jy/ Jx and Jz/ Jx- This VEV was found in work (see 
Appendix A, Eq. ( A.7 )). As for the two-particle form factors, the Z2-symmetry generated 
by (IllD enables one to predict their general form: 

(^■)( (To" I B±(A:i)B^(fc2)>i, = (-1)^' i^f (A:i, k^) ± F^ikiM) , (3.14) 
^^\cjI I B±(fci)B±(fc2)>i, = Fl{k^, k2) ± i-iy F^ik^M) , (3.15) 
(^■)( cTo^ I B±(fei)B±(fe2)):„ = (-1)^' Ff{k^, k2) ± F|(fci, fc2) , (3.16) 
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where ^'^^ some functions of the quasi-momentums ki and The invariance with 
respect to charge conjugation C ( ^.11 ) dictates that all other two-particle form factors vanish. 
Form factor ( 3.14| ) was calculated in Ref. [27|. It is possible to extend the result of that work 
and calculate form factors ( 3.15| ) and ( p.l6 ) as well. We collect the explicit expressions of F 
and iii Appendix A. 

To conclude this subsection, let us note that, in the current treatment of the infinite XYZ 
model as the thermodynamic limit of finite periodical chains with an even number of sites, 
we cannot construct excited states containing an odd number of the elementary excitations 



B 



...B 



vac)(i"J') G V(^') 



(m = 0,l...) . 



(3.17) 



Such states are deduced from finite chains with boundary conditions breaking the translation 
invariance. In particular, the linear subspace of V^-'^ spanned by the states ( 3.171 ) can be 
constructed from the thermodynamic limit of chains with an odd number of sites and 
subject to the so called twist boundary condition^: 



JV+l ) 



a 



N-l 
2 



O' N+1 

2 



The in-asymptotic states (|3.9| ), ( ^.17 ) form complete bases in V^-'^ for < 0. 



3.2 XYZ spin chain in the scahng Hmit 

For Jx = Jy, the gap in the spectrum of the XYZ chain vanishes, and its correlation 



length |29| 



Rr 



8 (Jx - Jz. 
Jx{Jx Jy 



1 

2-2)7 



1 



vr 



arccos{Jz/ Jx) 



(3.18) 



becomes infinite. In the limit Jx Jy, the correlation functions at large lattice separation 
(~ Rc) assume a certain scaling form which can be described by quantum field theory. If 
( Jj. — Jy) I Jx <C 1, it is natural to treat the XYZ model as the perturbation of the XXZ chain 
by the lattice operator 



The leading term in the expansion (|T]^) for this operator is given by the relevant field 
^2,0 + 0-2,0 of scaling dimension 2r/. Therefore, the scaling limit of the XYZ chain is 
described by the sine-Gordon quantum field theory |2£ 



Hsg = Hcauss - j dx (©2,0 + 0-2,o) (x) 



(3.19) 



Up to a numerical factor, which was obtained in Ref. [30|, the coupling constant fi coincides 
with the quantity M^~^^, where the combination 



M = {eRc 



(3.20) 



^This fact follows from the result of work 



13 



can be naturally identified with the soliton mass in the sine-Gordon model (3.19). 



The sine-Gordon theory admits a class of local soliton-creating operators characterized by 
two integers s,n S Z, where n gives the topological charge and sn/2 represents the Lorentz 



spin of the field. These operators can also be expressed in a form similar to ( 2.14 ), in which 



ip denotes the sine-Gordon field instead of the Gaussian field obeying the simple D'Alembert 



equation ( ^.S| ) (see e.g. Ref. |31| for details). They moreover coincide with the Gaussian 



fields ( 2.14| ) in the conformal limit /x ^ 0. Hence, with some abuse of notation, we will 



denote such soliton-creating operators in the sine-Gordon model by the same symbol C^s^rifj* 
To proceed further, one needs to draw a link between local lattice operators in the XYZ 
chain and local fields in the sine-Gordon model. Let us note at this point that local expansions 
of the type ( |1.2[ ) are based on dimensional analysis and do not necessarily imply the criticality 
of the original lattice system. Similar expansions are expected to be applicable to describe 
the near-critical behavior of lattice systems. Usually, the term with the smallest scaling 



dimension in (1.2) governs the universal scaling behavior of lattice correlators, whereas the 
next terms produce non-universal corrections. In particular, relations (|3.1| ), ( |3.2D obtained 
for the XXZ spin chain can be used to study the leading scaling behavior and first non- 
universal corrections of the XYZ correlation functions and form factors. In the XYZ case, 



the continuous fields which appear in (3T), ( p.2| ) should be understood as operators in the 



sine-Gordon model rather than their conformal limits. The numerical constants Cq, Ci and 
Ci remain, of course, the same as for the critical XXZ chain. 

Let us now discuss the relation between the Hilbert spaces of the XYZ and sine-Gordon 
models. In general, the theory ( 3.19| ) admits a discrete symmetry ip ^ (f + 27rj'/y^ {j G 



which is generated by the operator defined similarly to (2.19). For < 77 < 1, the above 
symmetry is spontaneously broken, so that the theory has an infinite number of ground 
states I Oj ) {j £ Z) characterized by the corresponding VEVs of the field (p: 

{Oj\p{x)\0,) _27ri 

(o,|o,) ^ ■ ^ ■ ^ 

The sine-Gordon model which governs the scaling behavior of the even sector of the XYZ 
chain is subject to the additional constraint = 1. This equation implies in particular that 
the field p is compactified, (f = (f + 47r/y^, and that, unlike the uncompactified case, there 
exist only two non-equivalent vacuum states | Oj ) with j = 0, 1. These states are naturally 
identified with the scaling limit of the two XYZ vacuums | vac)^-'^ 

To describe the scaling limit of XYZ excited states, one has to relate the Zamolodchi- 
kov-Faddeev operators of the lattice and continuous theories. Let us recall here that the 
sine-Gordon model admits a global continuous U{1) symmetry generated by the operator 

V„ = e^"^, where Q = ^ [ dxd^ip, (3.22) 



2vr 

which acts on the exponential fields as follows, 

®In Ref. the soliton-creating operators were denoted as O", where a — Syfrj/2. The parameter /3 

in coincides with y/rj. Note that there, the quantity 2a //3 was not assumed to be integer. 
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Notice that the Gaussian CFT also possesses such a global symmetry, contrary to the XXZ 
and XYZ lattice models. Nevertheless, the form of the expansions (|]3|), (|]3|) and ( ^) 
suggests that the lattice models are invariant with respect to the transformation which 
acts trivially on all local lattice fields. Such symmetry manifests itself in the existence of 
two subspaces V^-'^ {j = 0, 1) which can be treated as eigenspaces of the operator V^r: 

v^v(^') = (-1^ v(^') . 

The fundamental sine-Gordon kink-like excitations, the soliton A_ and the antisoliton A_|_, 
carry, respectively, negative and positive units of the topological charge Q ( 3.22| ): 

V«A±(0)V-i = e±^- A±(0) , (3.23) 

where the argument 9 denotes kink rapidity. The relation with the Zamolodchikov-Faddeev 
operators of the XYZ model was established in Ref. [32|. In our notations it can be summa- 



rized as follows: the operators A-|-(A;) defined in the lattice model as 

a£'°) = ^ ( ± bJ + bL''°) ) , Ai°'^) = 4 ( - B?') ± bL°'') ) , (3.24) 
\/2 \/2 

turn out, in the scaling limit, to be the Zamolodchikov-Faddeev operators A±{9) (|3.23| ) of 



the sine-Gordon model. Here, as well as in Eq. (3.11), we denote the restriction of the 



operators B-t and A± acting on the subspace V^-'^ as b||1 "''•'^ and a£ ■''"'^ Again, with 
some abuse of notation, we use the same symbol A± for the XYZ operators and for their 
scaling limits. Notice that the quasi-momentum k of the low-lying fundamental excitation 
becomes the usual particle momentum in the scaling limit: 

k £(k) 
hm - = M sinh(6') and hm = M coshie) , 
£^0 e e-*Q e 



where M is the soliton mass ( p.20 ). 



3.3 Scaling behavior of the form factors 



We are now in a position to calculate the constants appearing in expansions (3J.), ( |3.2D . To 
illustrate the procedure, let us consider first the scaling behavior of the lattice VEV ( p.l3| ). 
From relation ( |3.1[ ) we deduce that 



F~ 2Coe-2{Oifl) + ... . (3.25) 



Here (...) means an averaging with respect to the vacuum ( 3.21| ) with j = 0. To write the 



equation ( ^.25 ), we also use the fact that the VEV of the fields Oifl and 0-i,o are equal by 



virtue of charge conjugation symmetry. The VEV of the operator Osfi was found in ||3. 
We denote it as 

( Cs,0 ) = \/ Zsfi ■ 



On the other hand, the function F is given by the Baxter-Kelland formula ||2^ (see Eq.( A.?]) ). 
Its leading scaling behavior reads 



(4i?,)-2 +... . (3.26) 
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Comparing Eqs. ( p.25| ) and ( |3.26| ) , and using the relation R^^ = Me, one can deduce the 
constant Cq: 

C„= ^„ \ (3.27) 



2(1 - r,),/ZT^„ y 4 

A similar strategy can be applied to calculate the constants Ci and Cf in ( |3.1| ), p7 
The relations ( |3.24| ) allows one to express the leading scaling behavior of the functions that 
appear in Eqs. ( |3.14| )-( ^.16 ) through the form factors of the sine-Gordon fields: 

Ff ± i ~ 2 Co £2 ( O±i,o(0) I A+(ei) A_ (^2) >;„ + •••, (3.28) 

F| ± i F| ~ -Ci e'^^^ { O±i,2(0) I A_(0i)A_(^2) >;, + ■••, (3.29) 

and 

i^f (5t^(0)|A+(^i)A_(^2)L + ... , (3.30) 

Fir^'-Cfe^ (Oo,2(0) I A_(0i)A_(02) >;, + ... . (3.31) 



We use here an abbreviated notation similar to ( |3.12| ), except that the index specifying the 
vacuum states is omitted since it is always assumed to be j = 0. The two-particle form 
factors of the topologically neutral operators 0±ifi and have been known for a long time 
(see e.g. [25|). They read explicitly, 

( 0^1,0(0) I A+(^0A-(^2) = VZ^O .'V . t ' , ^ . , , (3.32) 

91+62 ei+92 X 

"'^ ^/r] G{-tTT) cosh [{9i -62 + n:)/{2£)) 

where ^ = ^.nd the function G{9) is a so-called minimal form factor. The form factors of 
the topologically charged operators Og^n have been proposed in In particular, for n = 2 
one has: 

( a,2(0) I A_(0i)A^(02) >i„ = e'i' 6^^+^ G(0i - 02) . (3.34) 

The explicit formulae for the minimal form factor G{9) and the field-strength renormalization 
Zg^n are recalled in Appendix B. 

We are now able to compare ( p.28| )-( |3.31| ) with expansions of the exact lattice form factors 
(|3.13| )-( ^T6|) given in Appendix A, and to relate the values of the constants Cq, Ci and Cf 
to the constants Zs^n'- 

C,= (3.35) 
Cl= , ^[^f. (3.36) 
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4 Correlation amplitudes. Comparison with numerical re- 
sults 



With the explicit expression ( |B.3| ) for the normahzation constants Zi^, Z\^2 and Zo,2; the 



relations (3.27 ), (p.35| ) and (3.36) lead to the following formulae for the correlation amplitudes 
A, A and A^ (p.SOl): 



A 



1 



2(1 -r/)^ 



2V^r( 



2-21)1 



exp 



sinh(r/t) 



t vsinh(t) cosh((l 



Ti)t) 



rj e 



-2t 



(4.1) 



A 



r/(l -r/) 



r( 



2-2r]' 



2-2i)l 



X exp 



dt 



cosh.{2r]t)e 



-2t 



1 



+ 



t \2sinh(r/t) sinh(t) cosh((l — r/)i) sinh(ryt) 



rf±l^-2. 



and 



2V^r(^: 



exp 



/■°° dt / sinh((2r? - l)t; 
Jo i Vsinh(?7t) cosh((l — 



r/)t) 



2r? - 1 ^-2t 



Note that these amplitudes obey the simple relation: 



A 

'AA. 



r2fi + 



TT 



^^\\^2- 



2-2i? 



) 



2r? 



(4.2) 



(4.3) 



(4.4) 



The correlation amplitude A was already obtained in [33|. Our computations also confirm 
the conjecture from concerning the amplitude A^. 

In ||3^, numerical values of the spin-spin equal-time correlation functions have been 
obtained for an open chain of 200 sites by the density-matrix renormalization-group tech- 
nique 11 . In Table I, we compare, for different values of the anisotropy parameter A, the 
numerical values that follow from ( f4.2D , (4.3) with those estimated in |^5| from the fitting of 
the numerical data. The corresponding plots are represented in Figure ^. 



5 Spin-spin correlation functions in the isotropic limit 

5.1 Marginal perturbations of the Wess-Zumino-Witten model 

As long as the parameter r] is not too close to unity, the first terms of the asymptotic 
expansions (2.47), ( 2.48[ ) provide a good approximation to the spin-spin correlation functions 
even for moderate space separations /. However, these expansions cannot be directly applied 
in the isotropic limit r] ^ \. Indeed, in this limit, the operator ©0,4 + 00,-4 in the effective 
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A 


7 AfTT M 1 A 
Ai\ u ivi n\ 


1/4 






0.7 




0.004720 


0.008(1) 





00893 


0.6 


0.0048(14) 


0.006643 


0.0133(1) 





01314 


0.5 


0.0076(9) 


0.008656 


0.0184(4) 





01795 


0.4 


0.0099(7) 


0.010696 


0.0235(2) 





02332 


0.3 


0.0122(4) 


0.012717 


0.02921(3) 





02924 


0.2 


0.0144(2) 


0.014687 


0.03556(3) 





03574 


0.1 


0.0164(2) 


0.016583 


0.0425(2) 





04285 


0.0 


0.0182(2) 


0.018386 


0.0501(5) 





05066 


-0.1 


0.01995(7) 


0.020082 


0.0588(3) 





05929 


-0.2 


0.02154(5) 


0.021657 


0.0683(6) 





06891 


-0.3 


0.02296(4) 


0.023098 


0.0791(8) 





07978 


-0.4 


0.02420(4) 


0.024392 


0.0918(9) 





09231 


-0.5 


0.02525(6) 


0.025522 


0.1063(9) 





10713 


-0.6 


0.0261(2) 


0.026464 


0.1236(5) 





12539 


-0.7 


0.0267(3) 


0.027182 


0.145(1) 





14930 


-0.8 


0.0271(7) 


0.027608 


0.171(5) 





18414 


-0.9 


0.027(2) 


0.027570 


0.20(1) 





24844 



Table 1: Correlation amplitudes. A and Az are computed from expressions (^), (|4.3|) , 
whereas A^^^ and A^^^ were obtained from the fitting of the numerical data in |^5| . The 
number in parentheses indicate the error bar on the last quoted digits. 



Hamiltonian ( 2.40 ) becomes marginal and induces logarithmic corrections to the leading 
power-law asymptotic behavior. 

The suitable way to explore the — > 1 limit is based on the low-energy effective theory 
defined as a perturbation of the Gaussian model with r/ = 1. As is well known (see e.g. 
|17|), the Gaussian model coincides in this case with the SU(2) level one Wess-Zumino- 
Witten (WZW) theory. In particular, the WZW holomorphic currents are identified with 
the following primary operators of the Gaussian CFT, 

rR = \{dt-d,)^, J% = 0±i,^2, (5.1) 

Jl = ^(dt + d,)^, J^ = 0±i,±2, (5.2) 

whereas the matrix of the fundamental WZW field is bosonized as 
00,2 iO-ifi 

The low-energy effective Hamiltonian can be expressed as a marginal current-current per- 
turbation of the WZW Hamiltonian |37|, 

Hxxz = Hwzw + I ^{9\\JRJl + Yi-^BJL +JrJI) + ■■■}■ (5-3) 
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Figure 1: The correlation amplitudes A and as functions of the anisotropy parameter A. 
The points with the error bars represent the numerical data obtained in (see Table 
and the continuous line follows from Eqs. ([4.2| ) and (|4.3| ) . 



In this expression, the coupling constants ^ly and g± should be understood as "running" 
ones, i.e. depending on the renormalization scale r which has the dimension of length. The 
corresponding Renormalization Group (RG) flow is known as the Kosterlitz-Thouless flow. 
For our purpose, we need only to consider the domain 

\9±\^9\\, (5.4) 

in which all RG trajectories flow toward the line g± = the infrared-stable fixed points 
associated with the Gaussian CFT. These trajectories are characterized by the limiting values 
of the running coupling g^^, 

e = J lim 5^11 (r) , (5.5) 

and the parameter e is simply related with the parameter ry of the Gaussian model, 

e = l-7?. (5.6) 

The RG flow of the running coupling constants is defined by a system of ordinary differ- 
ential equations, 

^ ^ gj ^d9± ^ 9\\g± 

dr f\\{g\\,g±) ' dr f±{g\\,9±) 

Perturbatively, the functions f\\^±ig\\, g±) = 1 + 0{g) admit loop expansions as power series 
in 5|| and g±, and their precise form depends on the choice of a renormalization scheme. We 
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use here the scheme mtroduced by Al.B. Zamolodchikov [30,3^, who showed that, under a 
suitable difFeomorphism in g\\ and g^, the functions /y and /_|_ can be taken equal to each 
other and to the quantitj^ 

/|| = /± = 1 - I . (5.8) 

With this particular choice of the /3-function, it is possible to integrate the RG flow equations 
exactly. To do this, let us first note that the system of differential equations (|5.7D , (|5.8|) 
admits a first integral, the numerical value of which is determined by means of the condition 



(ED: 

gl-gl = {2ef. (5.9) 



Then the equations (|5.7|) are solved as 

51, =26^^, g^=4e-^, (5.10) 

where q = q{r) is the solution of 

1 1 , , > / \ ^ 



q2-.-2 (i-q)=et^^j . (5.11) 

As well as e, the dimensional parameter rg is a RG invariant. It is of the same order as the 
lattice spacing e, and is supposed to have a regular loop expansion of the form 

— = exp(co + cie + C2e2 + ... ) . (5.12) 

ro 

It should be noted that the even coefficients cq, C2, ... in ( |5.12| ) are essentially ambiguous 
and can be chosen as one wants. A variation of these coefficients corresponds to a smooth 
redefinition of the coupling constants which does not affect the /3-function. On the contrary, 
the odd constants C2k+i are unambiguous and precisely specified, once the form of the RG 
equations is fixed. It is possible to show [^,^| that the odd constants vanish in Zamolod- 



chikov's scheme: C2k+i = (/c = 0, 1 . . . ). Therefore, once the coefficients C2k in ( 5.12| ) are 
chosen, the renormalization scheme is completely specified. 

As already mentioned, the perturbation by the marginal operators produces logarithmic 
corrections to the scale-invariant correlation functions. Hence, the conformal normalization 
condition imposed on the field Os,n turns out to be singular for e = 0, and we should define 
renormalized fields 0^^^\ which are rescaled version of the "bare" exponential operators: 

0(7)(a;,t;r) = Z~l{r) Os,n{x,t) ■ (5.13) 

Notice that, in writing (|5.13| ), we assume that there is no resonance mixing of the operator 
Os,n with other fields, so that it is renormalized as a singlet. In particular, one can easily 
check that this is indeed the case when 

|n|<2 + 2|s|. (5.14) 



""The analysis of Refs. Q concerns the ultraviolet-stable domain of the running coupling constants, 
which in the current notations is defined as — > \q_L\. Due to the perturbative nature of the RG flow 



equations, the /3-function for the infrared-stable domain ( 5.4) can be obtained through the formal substitution 
9\\ — * ~9\\ in the original Zamolodchikov equations. 
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The renormalized fields ( 5.13| ) are no longer singular at e = 0, but depend on the auxihary 
RG scale r. To specify them completely, we have to impose some non-singular normalization 
condition. The conventional condition, which is usually imposed on Green's function for a 
space- like interval — < 0, is 



(TO(X^(x,t;r)0(X^t(o,0;r)>|^ 



t — X 
t + x 



(5.15) 



Eqs. ( ^.13| ) and ( ^.15| ) imply that the correlators of the "bare" exponential fields should take 
the form: 



{TOs,n{x,t) Ol,„(0,0); 



t — X 

t + X 



(5.16) 



Then, using (3T), (|3.2D , one can express the spin-spin correlation functions through the 
renormalization factors Zg^n- For example, for the time-ordered correlation function of cj^, 
one has: 



^(^ + ^)-7^(e0;^) ^ + .... (5.17) 



Here we use notations ( p. 49 ), and the function IZ is related to the following causal Green's 
function as, 



( T Oi,2(x, t) O_i,2(0, 0) > = Zi,2(r) 7e(r 



r=\/ x'^—t'-^ 



(5.18) 



The first terms of the perturbative expansion for the scalar factor Zg^ni'i') in ( ^-13 ) can 
be deduced from the results of work [pyl: 



^.,n(r) = ( -) ^ (<7i)^"~^ ^ 



X e 



u\g\\+u2gr. 



(l + 9l{vi-V2g\0+O{g^)), (5.19) 



where 



and ds^n is given by ( p. 15 ) with r/ = 1 — e. The coefficients ui, U2, vi and t>2 in these 
equations are listed in Appendix C. Notice that to derive ( 5.19| ) one should assume that the 
field Og^n is mutually local with respect to the density of the effective Hamiltonian ( |2.40| ) . 
This assumption implies that s S Z, but does not impose any restriction on n in addition to 

(ill). 



As follows from the Callan-Symanzik equation, the function IZ in ( 5.1 j ) admits a per- 
turbative expansion in terms of the running coupling constants. Explicitly, one can obtain 



7^ 



T 1 ^1 



+ {c-\)gl + cgl + 0{g')Y 



(5.20) 
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The constant c appearing in Eq. ( 5.20| ) is related to cq from ( 5.12| ) as 

CO = c + 7s + ^ ln(27r) , (5.21) 

where = 0.5772... is the Euler constant. Combining relation ( 5.17] ) with ( ^.191 ) and 
(l5.20| ), one can deduce the RG improved expansion of the lattice correlator which is 
applicable for e <C 1. 

We can similarly derive an expansion for the correlation function of . The relation 

dt^ = — dAOo,n , (5.22) 

2tT^ ITT n=0 

which follows from the definition ( |2.14| ), is useful to perform this computation. Note also 
that the operators 

Co,2 + ^0,-2 C>o,2 — C>0,-2 



V2 ' 
renormalize as singlets: 

Indeed, since 

c o± c = ±o± , 



V2i 



(5.23) 



invariance with respect to the charge conjugation prevents resonance mixing of 0+ and O- 
Now, using equations ( 15^221 ), (K^) and (pD, one obtains 



(raf+,.(t)a|(0)) ~ i^l^+^l_f^lZo,n{e^/hl^) 



n=0 



(5.24) 



We collect in Appendix D the RG improved expansions of the different-time two-point 
correlation functions ( ^.17 ) and ( ^.24)) . 

5.2 Equal-time correlation functions for the XXX spin chain. Comparison 
with numerical results 



Using expansions ( D.l ), ( D.2 ), it is easy to perform the isotropic limit. Setting e = 
and g± = = g, one obtains the following large I expansion for the equal-time spin-spin 
correlation functions^ 



1 



7r3 

21 

+ 



1 + 



9 + 



128 16 



+ 



7c 7c2 c3 13C(3)\ 



1024 ' 256 64 16 32 
(-1)' f 1 , 5 , / 3x 5\ c(c + 2) 



IT- 



2/2 



1 + f + (C+T ) ^ + 



) + 0(/ 



(5.25) 



^The coefficient ^2/7r-^ in Eq. (5.25) was originaUy obtained in Ref. pof . 
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Here g = g{l) is a solution of the equation 



2V27re 



7E+C 



(5.26) 



which corresponds to the hmit e — > of Eqs. ( 5.10| ) and ( 5.11 ). Let us stress that, if the 
perturbation series in ( ^.25 ) can be summed, then the correlation function should not depend 
on the auxiliary parameter c (5.12), ( 5.21| ). This, however, is not true if we truncate the 
perturbative series at some finite order. Thus, when fitting the numerical data with ( ^.25| ), 
we may treat c as an optimization parameter, allowing us to minimize the remainder of the 
series, or at least to develop some feeling concerning the effects of this remainder. 

The correlation function ( ^.25| ) has been computed numerically in |43] for 1 ^ / ^ 30. 
The authors used the density-matrix algorithm to study the large-distance decay of the 
correlation function for XXX spin chains with 14 < < 70 sites. To extract the values 
of the correlation function in the infinite chain case, they adopted the phenomenological 
scaling relation of Kaplan et al. |^l| (see also |^2|). The relative error of the interpolation 
procedure was estimated to be of order 1% for the largest I values. In Table ||, we compare 
those numerical data to the results obtained from ( 5.25| ) in the cases c = — 1 and c = — 2. 
The corresponding plots (numerical data against RG result for c = — 1) are given in Figure 
^. It appears that the numerical data are consistent with our prediction within the given 
errors. 



I 


NUM 


RG (c = -1) 


RG (c = -2) 


I 


NUM 


RG(c= -1) 


RG (c = -2) 


1 


0.1477 


0.1616 


*** 


16 





1628 


0.1624 


0.1630 


2 


0.1214 


0.1213 


0.1918 


17 





1676 


0.1666 


0.1671 


3 


0.1510 


0.1509 


0.1583 


18 





1646 


0.1642 


0.1647 


4 


0.1384 


0.1381 


0.1424 


19 





1689 


0.1679 


0.1683 


5 


0.1541 


0.1541 


0.1566 


20 





1661 


0.1657 


0.1661 


6 


0.1463 


0.1462 


0.1482 


21 





1700 


0.1690 


0.1694 


7 


0.1567 


0.1571 


0.1586 


22 





1674 


0.1670 


0.1674 


8 


0.1513 


0.1514 


0.1527 


23 





1712 


0.1700 


0.1704 


9 


0.1596 


0.1596 


0.1607 


24 





1687 


0.1682 


0.1686 


10 


0.1550 


0.1552 


0.1561 


25 





1723 


0.1710 


0.1713 


11 


0.1620 


0.1618 


0.1626 


26 





1699 


0.1693 


0.1697 


12 


0.1581 


0.1581 


0.1588 


27 





1734 


0.1719 


0.1722 


13 


0.1641 


0.1636 


0.1643 


28 





1710 


0.1703 


0.1706 


14 


0.1606 


0.1604 


0.1611 


29 





1746 


0.1727 


0.1730 


15 


0.1659 


0.1652 


0.1658 


30 





1722 


0.1712 


0.1716 



Table 2: Numerical values of the correlation function ^ ( crfj^j crj ) of the XXX spin chain 
according to the distance I. The column "NUM" has been obtained in |^^, whereas the 
columns "RG (c = —1)" and "RG (c = —2)" follow from ( 5.25| ) with the corresponding 
values of the free parameter c. 
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0.18 
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0.14 
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W(l) (numerical) 
W(l) (RG) 




Figure 2: The correlation function W{1) = j i'^f+j ^ J ) of XXX spin chain according to 
the distance /. The dots represent the numerical data obtained in [43|, and the continuous 
line connects odd and even terms obtained from Eq. (|]25|) with c = -1. 



5.3 Erratum of 112 



The spin-spin correlation function in the limit e — > was previously studied in Section 5 of 
Ref. |jl2|. The analysis performed in that paper, along with the numerical results obtained 
in IQ (see Fig. 2 of |12|), strongly suggested the existence of an additional staggered term 
of the form 



oc 



-1)' + 1 



in the large distance asymptotic expansion of the correlation function ( af^j '^j )■ 
argued in [12| that such a term occurs because of the presence of a correction of the type 
(— 1)^ dxO±ifi in expansions (|3.1| ). However, the RG computation from appears to be 
erroneous. Indeed, Eqs. (5.6) and (5.7) from |12| have to be replaced respectively by our 
equations (D.l) and ( ^.251) . Therefore, contrary to what was claimed in the numerical 
data are consistent (within the numerical errors) with (|3.l| ). 



6 Conclusion and further remarks 



The purpose of this article is the quantitative study of the long-distance behavior of spin-spin 
correlation functions for the XXZ Heisenberg chain in the critical regime. Our main result 
here is the determination of analytical expressions for the correlation amplitudes involved in 
the corresponding asymptotic expansions. To obtain these values, we considered quantum 
field theory which describes the scaling limit of the lattice model, and compared, in this 
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limit, the respective normalizations of the lattice operators and of the corresponding local 
fields. This comparison was achieved by considering known exact matrix elements (form 
factors). 



We would like to conclude the article with the following remarks. 

• The method used in this work can be applied to higher order terms in expansion ( |2.35| ) 
or (|2.36D . For example, we were able to compute (up to sign factors) all constants Cm from 
( 2.36| ) for odd integers m = 2p +1: 



2p+l) 



r?(l -r?) 
X exp 



r( 



2-2ri' 



20Fr 



1 



°° dt 



n 

cosh(2r/t)e-2(2p+i)t 



sm" — J cot 



1 



vr(2j - 1) 
2-2r] 



t \2sinh(?7t) sinh(t) cosh((l — r7)t) 
2p + l 



+ 



sinh(ryt) 



(2p + l) X ^2t 
r] H I e 



(6.1) 



One can also study expansions of other local lattice operators in terms of the scaling 



fields. For example, we have calculated the constant Cq 
of the lattice operators 



is) 



in the leading term of the expansion 



eye 2 



for which we obtained (p = 0, 1 . . . ): 



a 



(2p) 



^ \2-' 



2-2ri' 



G 



(2p+l) 




2^/^r(2^) 



2r]p^ 



2^T{^. 



7rP(l - rifP^^P 



2r,(p+i)2 



1 



r2/ '?(2j-i) 

^ V 2-2»7 



(6.2) 



(6.3) 



2-2J7 

X exp 



27rP(l 
°° dt 



^)2p2+p+l 







sinh(?7t) 



y2( m 



!Z e-2* 

t V 2 sinh(t) cosh((l - r/)t) 2 



(6.4) 



• Eventually, one can wonder if it is possible to confirm our predictions from existing 
integral representations of lattice correlators. Up to now, although explicit expressions for 
the equal-time spin-spin correlation functions at finite lattice distances are known [Q, Q , 
their long-distance behavior was studied only for the so-called "free fermion point", A = 0. 
In this case the XXZ spin chain can be mapped onto two non-interacting critical Ising models 



and the long-distance asymptotics are readily derived from results of works p5|-|4^ . 

• The approach of this work is actually quite general for lattice solvable models at criti- 
cality: from a knowledge of particular form factors of a lattice theory and of its quantum field 
theory counterpart at the scaling limit, it is possible to predict the amplitudes which govern 
the large distance behavior of lattice correlation functions. It would indeed be interesting to 
obtain effective results for other critical exactly solvable lattice models. 
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A Two-particle form factors in the XYZ model 

In this appendix, we collect explicit expressions of two-particle form factors of local spin 
operators in the XYZ model. Following Baxter fs], we use the parameterization of the 
coupling constants Jx > Jy > \ Jz\ the Hamiltonian ( |3.3| ) in terms of < < 1 and of the 
elliptic nome < p < 1: 



(Mn)fM_ _ ^i(^Ki(o) ' 

vre V^4m(7?) i?4(0)T?4(r/); ' 



, _i-v f^'Av) ^'M\ 

Here 'di{u) = ■di{u,p) denote the elliptic theta- functions 

oo 

??i(n,p) = 2pi/^ sin(7r'u) _ 2/" cos(27rn) + /") , (A.4) 

n=l 

oo 

i?4(n,p) = _ 2p2"-icos(27rn) +/(2n-i)) , (A.5) 



n=l 



and the prime in Eqs. ( A.l| )-( ]Al3 ) means a derivative: 'd\ = du'&i- We shall also use the 



other conventional theta-functions 

Mu) =Mu + 1/2) , Mu) = + 1/2) , 

and the notation, 

C = . (A.6) 

1 — r] 



With this parameterization, the VEV of ( p. 13 ) is given by the Baxter-Kelland for- 



mula [g8|: 
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In order to describe the two-particle form factors, one needs to know the exphcit form 
of the dispersion relation ( 3.10| ). For this purpose, it is convenient to parameterize the 
quasi-momentum k by means of the so-called rapidity variable 0: 

e^'=(^) = ^^(^-i'P;^) . (A.8) 
As a function of 6, the excitation energy explicitly reads [ p^ , 



Equations ( [A.8| ) and ( |A.9| ) define the dispersion relation £ = £{k) in parametric form. 

The two-particle form factors of the spin operators can be computed by means of the q- 
vertex operator approach, for which progress has been made recently in the XYZ case j26,27|. 
In [^] , the two-particle form factors of the operator were obtained]^: 



Fq Giei-92,p) ^4(0, M'-^^p "^) 

r,^\ „0 ^ 01-02+iTT 



(A.IO) 



' Mit-lp'^)M§,-hp-^) M'-^^p"^)' 

For —2tt< Qm{6) < the meromorphic function G reads 

.2/ 



^ ^^Jl±f (e+^.)2 r y 1 sm\6nie + i7r)/2) sinh(7r<5(e + l)n/2) 
' I sinh(7r(5n) sinh(7r(5^n/2) cosh(7r5n/2) 



(A.12) 



and it is defined through an analytic continuation outside this domain. The parameter 6 in 
( A.12| ) is related to the elliptic nome as p = e , and the constant Fq in ( [A.IO ), ( A.ll| ) 

is given by 

^0 — :r- ; TT ■ 

< e[{o,p2) 

Notice that, in writing the form factors, we always assume the conventional normalization 
of vacuum states, ( vac | vac ) = 1, and of in-asymptotic states: 

n 

in( B^,M ■ ■ ■ ^^[(k'l) I B.i(A;i) . . . )i, = (2vr)" J] ^-.-^ ^(^^^ " ' (^.13) 

i=i 

where kj = k{6j) and kj = k{9j). 



^Note that the regime considered in is the so-called principal one { — Jz > Jx > \ Jy\)- To apply the 
results obtained there to the case with Jx > Jy > \Jz\, one has to replace af , af and erf from respectively 
by af , (— and (— l)'(Tf , which corresponds to a similarity transformation of the Hamiltonian (3.3). 
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Using the method proposed in |27|, one can also compute the two-particle form factors 
of the other spin fields, and a^: 



py^_ FoGi9,-e,,p) ^3(0, P^) M'-^^P-^) ,A 14) 



j,y^_ FoGie,-e,,p) ^3(0, P^) M'-^^p^) .A 15) 

' Mit-k^p'^)MA-k,p-^) M'-^^p"^)' 



and 



M^-Ip'^)M^.-\,p'^) M'-^^p-^) 



(A.16) 



p.^^ FoG{e,-9,,p) ^2(0, M'-^^p"^) .A17) 

To compare the lattice and the sine-Gordon two-particle form factors, one should take 
the limit p ^ 0. Notice that 

where the correlation length is defined as in ( P^ ). The following relat ion between the 
function G ( |A.12| ) and the sine- Gordon minimal form factor 

lim G{e , p) - '^^^^ 



P^o ' G{-iTr) ' 
is useful to proceed with this limit. 



B Form factors of topologically charged operators in the sine- 
Gordon model 



In this appendix, we recall the expressions obtained in ||31| concerning the form factors of 
the topologically charged (or soliton-creating) operators Os,n in the sine-Gordon model. The 
simplest non-vanishing form factor of Os^n is given by the formula: 

n 

(a,n(O)|A_(0i)...A_(^„)>.^ = V^e^ n n G{e„,-e,). (B.l) 

m=l rn<j 

Here, the minimal form factor G has a form: 

^/.x ^ w.,, N f /""^ c^t sinh2(t(l-i0/7r)) sinh((e-l)t) 1 , 

^(6* = iCi sinh 61/2 exp<^ / \\ , ' ' , . \ ■ B.2 

^ ^ ^ ' ' \Jq t sinh(2t cosh t sinh^t J ^ ' 
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The explicit expression of the normahzation constant Zg^n^ which has been conjectured in |31], 
is the following: 



(--2 \ 2 / ? (--2 \ T 



^Mr(| + |) 



3 I i\ -1 2d3,n 



X exp 



ICl) V 16 

°° dt r cosh(2^st) e-(i+«)"* - 1 







t 



+ 



n 



2 sinh(^t) sinh((l + i)t) cosh(t) 2 sinh(tO 
In the previous formulae we use the notations, 

sinh2(f/2) sinh((^ - \)t) 



2 d.q.r}. G 



Ci = G(— ivr) = exp 



t sinh(2t) cosh(t) sinh(^t) 



C2 = exp <^ 4 



dtsinh2(t/2) sinh((^ - l)t) 
~t sinh(2t) sinh(^t) 



and ^ is given by (|A.6|) . 



(B.3) 

(B.4) 
(B.5) 



C Numerical coefficients for equation ( [b.i9|) 



We collect in this appendix the explicit expressions of the coefficients ui, U2, vi and V2 which 
occur in the expansion ( [5.19| ). 



U2 



3\ s{s-\) 
16 V °\2) 2/ 4 
(n^ - 4^2) (n^ + 4s2 - 



3072 



n 



^11 1,0^ n(n?—4:)./n 



+ 



3n^ 



192 



n 



s{s + 2) lln^ c C2 (n^ — 
192 768" ^ 24 32" 



As' 



where 



Ts{z) = ipiz + s)+ ipi-z + s) + 2'yE + 2c , 
T',iz)=d,Tsiz), T':{z)=dlTs{z), 



(C.l) 
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and ip{z) = dzlogT{z). The constants c and C2 are the same as in Eqs. ( p.l2 ), (5.21). The 
coefficients vi and (using the expressions ui^2 from ( |C.1|) ) V2 are 



Vl 



V2 



128 ^2/^ 
3ra^ + 4g(2 - 5g) 
64 



64 



V2 

.2 



^ n\ + 4s(10 - 17s) ui 

2 J 128 T ' 



(n2-4s2)(8-3n2) ^„/n 
3072 ' V 2 



n(n^ - 4s^ 

128 
- 4s^ 



n(n2 + 45^ - 8s 
256 



T 



192 ^2 



n 



+ 4s{s -2) 2/n\ - 8(s2 - s + 1) 



128 



128 



^^^^(2c2-14C(3)-3)-^^^ + ^ + ^ + ^-^ 
128 ^ ' ' 64 8 2 2 8 



(C.2) 



D Spin-spin correlation functions for e <C 1 

In this appendix, we give RG improved expansions for the different-time spin-spin correlation 
functions which were discussed in Section ^: 



( rar+,.(t)aj(0) ) ~ J- — 



X g«?9||+«2Si] 

9 



{i + 5iK-^f5i,) + o(/)}-i^^ 



.2 )x e-(-^B+f logCS'^))^^ 



X 



xexp^ J + - + - 5 -(tttt---- 



2 ' 2-^-^ ' V96 
1 //+ l- 



+ 



and 



I- . , , ,0 



hi- (1 - f ) 



2 /, 



2c^ + c-i)^ + 0(/))+f(l + f2c-i 



2j5||+0(r 



(D.l) 



+ .... (D.2) 
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In these expressions, the constants are given by 



and 



W2 



Ml 



3 c 
8 ~ 2 ' 

~32 8 ~ T 



wi = c, 

1 c(l + c) 



8 2 
3 c 
8^2' 

5 5c 3c2 
~32 



U2 



W3 

Un - 



1 C(3) C2 

64 48 8 ' 

— --C(3) + ^-^, 
128 96 ^ ^ 6 8 



1 



13 



32 24 

64 48 8 8' 

11 71 , , c 5c2 2 o C2 

+ — C 3) + + — 

128 96 ^ ^ 4 4 3 8 



The running couphngs g^\, g± in ( p.l] ), (p. 21) are defined by equations ( 5.1C| ) and ( 5.11 ) where 

t 



r 



x/27r/+/_ e 



and 



l± 



L 



Setting 
one obtains (5.25). 



I (equal time) and g± = g\\ = g, e = (isotropic hmit) in ( p.l[ ), ( p.2| ) 
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